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Abstract 

We consider the effect on matter-enhanced neutrino flavor trans- 
formation of a randomly fluctuating, delta-correlated matter density. 
The fluctuations will produce a distribution of neutrino survival prob- 
abilities. We find the mean and variance of the distribution for the 
case of solar neutrinos, and discuss the possibility of placing a limit 
on solar density fluctuations using neutrino data. 

1 Introduction 



Matter-enhanced neutrino oscillations, especially in connection to the solar 
neutrino problem have been extensively studied in the recent years. More 
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recently some interest has developed in the problem of neutrino flavor trans- 
formations via the Mikheyev-Smirnov-Wolfenstein (MSW) effect in a ran- 
domly fluctuating matter density A general approach to neutrino oscillations 
in such inhomogeneous matter was developed in Ref. ||. A Boltzmann-like 
collision integral with blocking factors, describing the decoherence of neu- 
trinos in matter, was given in Ref. 0. Matter fluctuations which are not 
random, but harmonic [|], [|] or changing stepwise || were also considered. 

Redfield equations for a neutrino traveling in a region with delta-correlated 
Gaussian noise were recently developed in Ref. and applied to two-neutrino 
flavor transformations in the post-core bounce supernova environment in 
Ref. In parallel to these papers, an analytical procedure to calculate 
the survival probability was described in Ref. ||, and further implications of 
solar matter density random noise upon resonant neutrino conversion were 
studied in Ref. JTO] . 

The aim of this paper is to expand the analysis of Ref. J7| to investigate 
the mean and variance of the distribution of neutrinos when a randomly- 
fluctuating, delta-correlated electron density is present in the sun. A general 
treatment of fluctuations is presented in Section |2|. Mean survival proba- 
bilities and the variances of the survival probability distribution are given 
in Sections § and f|, respectively. In Section 5, we discuss the results and 
present conclusions. 
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2 General Treatment of Fluctuation 



We will be concerned with systems whose evolution is described by the 
Schrodinger-like equation 

i^(t)=H(t)m, (i) 

where ip is represented as a column vector. Following the method of Loreti 
and Balantekin |7j], we define the density matrix 

p = <0®^t (2) 

and divide H into two parts, one with known time dependence and one which 
fluctuates with time: 

H(t) = H (t)+B(t)M, (3) 

where B(t) is a c- number and the operator M does not depend on time. 
We assume that the fluctuation B(t) obeys 

(s(*o> = o 

(B(*i)fl(t 2 )) = a 2 f 12 
(B(t x )B(t 2 )B(t 3 )) = 
{B^BfoWhWu)) = a\f 12 f u + Aa/M + / 14 / 23 ) (4) 

where /jj = — gives the correlation between fluctuations in different 
places. Throughout the current paper, we will consider the case of delta- 
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correlated (white) noise: 

f(x) = 2t5(x), (5) 

with the correlation time r as a parameter. This is equivalent to the state- 
ment that the probability of a given B(t') is proportional to 

f exp{[B(t f )] 2 /2ra 2 } dt'. (6) 

J 

Further, the results from delta correlations will be approximately the 
same as those from step-function correlations of step length r as long as we 
have the constraint 

T<< ( C{t) + Jt hga{t) ) ' (?) 
where C(t) is the largest element of the matrix commutator [Ho, M\. Then we 
have the following result for the evolution of fluctuation-averaged p (cf. Eq. (16) 
of Ref. |0|): 

!<p(t)> = -a 2 r[M, [M, (p(t))]\ ~ W), (p(*)>]- (8) 

MSW conversion between neutrino flavors obeys the equation 

.d_( v e \ _ Am 2 / ((t) - cos26 sin2# \f v e \ , . 

l dt \ v x ) ~ AE \ ^ sin2^ -(C(*) -cos20) )\v x ) { ' 

where 

m = ^f^«w> do 

6 is the vacuum neutrino mixing angle, Am 2 is the difference in the squared 
masses of the two neutrino species, E is the neutrino energy, and N e is the 
electron number density. 



In order to study the influence of density fluctuations on the conversion 
rate, one can add periodic matter density perturbations to the average den- 
sity § 

JV e (r) = JVe[l + esin(Ar)], (11) 

where the wavenumber k is fixed. Such an additional term can induce ad- 
ditional MSW level-crossings giving rise to interference between them. To 
elucidate this behavior, one can utilize logarithmic perturbation theory, valid 
for small vacuum mixing angles. The application of logarithmic perturbation 
theory to the neutrino mixing problem is sketched in the Appendix. If there 
is more than one MSW resonance point, one can calculate the integral in 
Eq. (|5oD to obtain the electron neutrino survival probability, P e , as 



P, 



cxp 



Ti Am 2 sin 2 9 
4£ 



*a Am 2 
o 2E 



[C(t') -cos29)dt' 



;i2) 



If two turning points are close enough, one can utilize the uniform Airy 
approximation to obtain 



exp 



Am 2 sin 2 29 



' ix Am 2 r f b 



[C(t) -cos20] dt 



(13) 



4E Jt a 

which is the wavenumber dependence observed in Ref. [|J. 

If the density fluctuations cannot be parameterized with a single wavenum- 
ber as given above, but with a distribution of wavenumbers, one can then 
write 

N e (r) =Njr)[l + J dk e(Jfe) sin(Jfer)]. (14) 



The periodic fluctuation with a given wavenumber ko can be recovered by 
setting e(k) = eS(k — k ). The density fluctuations given in Eqs. (|5|) can then 
be considered as resulting from stochastically distributed e(/c)'s. 

3 Calculation of the Mean Survival Probabil- 
ity 

We will assume that the electron density N e fluctuates around the value N e , 
given by the Bahcall-Pinsonneault Standard Solar Model (SSM) including 
helium diffusion: 

N e = (1 + (3)W e (15) 

where (3 fluctuates and obeys constraints similar to those of Eq. (||) with 
delta correlations. 

We may then use the formalism of Section |2| with 

H = a z A(t) + a x B, M = cx 2 , (16) 

where we have defined 

A(t) = -^(C(t)-coB20), B = — sin20. (17) 
Then Rl we have the result 
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(18) 



k = 2(p 2 )r, (19) 
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and 



r = 2(z/> e )-l 
x = 2Re(^ e ) 

y = 2Im(z/> e ). (20) 



The condition of Eq. (|7|) now becomes 



,2\ - 1 



r«(sin2^j (21) 

where we have assumed that logiV e varies slowly with time. The right-hand 
side is the oscillation length of neutrinos at the resonance, divided by 4tt. This 



condition is similar to that of Ref. JTOl] that r <C A m , the oscillation length 



of neutrinos at any given point. For defmiteness, we will take the correlation 
length t to be 10 km. Then the constraint ([7]) becomes (assuming that the 
logarithmic derivative is small, which is accurate for the sun) 

sin 29 Am 2 /E < 3.95 x 10" 5 eV 2 /MeV. (22) 

We will present some results for which this condition does not hold, that is, 
for which the delta correlations we assume would give a different result from 
step-function correlations. The advantage of a constant correlation length is 
that the parameter k is constant in Am 2 / E, so that we can more meaningfully 
compare the effects of fluctuations with different MSW parameters than if 
we let r vary to satisfy Eq. (0). In any case, since r and (/3) rms enter only 
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through k, one can get identical results to ours by making r smaller and 
increasing (f3) 

rms • 

We solved the system (|X8|) numerically; results are presented in Figures 
1-3 for sin 2 29 equal to 0.01, 0.1, and 0.7, respectively, and a several values 
of (/9)rms- 0.01 and 0.7 correspond roughly to the small- and large-mixing so- 



lutions [0 to the solar neutrino problem. We used the standard solar model 
of Bahcall and Pinsonneault, including helium diffusion [|lTJ , but assumed 
all neutrinos were produced at the center of the sun. Usually, fluctuations 
suppress the v e — > transition. For the large-angle solution (sin 2 26 l = 0.7), 
the effect on the mean survival probability is noticeable only in the adiabatic 
region, even in the physically unreasonable case (/5) rms = 0.5. For smaller 
angles, there is some effect in the non-adiabatic region, but the greatest effect 
is still in the adiabatic region. 

In the adiabatic region, for survival probabilities greater than 1/2, fluc- 
tuations tend to enhance the transition. This has been explained |F] by 
appealing to the flatness of the Bahcall-Pinsonnealt density profile near the 
center. Neutrinos in this region do not go through the resonance point, but 
do go through part of the resonance region. Travel through that region has 
a significant effect since the density profile is flat just after production at the 
center. Fluctuations can therefore have an effect. To confirm that the en- 
hancement depends on the flatness of the density profile, we have also solved 
the system of equations (|T8| ) for an exponential density profile with central 
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density 6.25 x 10 25 cm -3 , and scale height 6.21 x 10 9 cm. These values were 
chosen to give the same density as the Bahcall-Pinsonneault model at the 
center and edge of the sun. At r = R Q , we cut the density to zero for both 
profiles. The exponential results are compared to the Bahcall-Pinsonneault 
results in Figures 4-6. For fluctuations above about 2%, the exponential pro- 
file gives noticeably weaker enhancement of the transition probability than 
the Bahcall-Pinsonneault profile, and changes the shape of the curve else- 
where as well. This result indicates that large fluctuations can induce appre- 
ciable non-adiabatic effects in the adiabatic region, since the initial density 
is the same in both cases, and both have the same numerical value of final 
density before the truncation to zero. (In any case, the step should have no 
effect, since the density at the edge of the sun is much less than the resonant 
density.) 

4 Higher Moments 

The formalism developed in Ref. j7| may be used to calculate higher moments 
of the distribution of survival probabilities. We will present results for the 
variance of the distribution, 

a 2 = (P 2 ) - <P e > 2 . 

In principle, the formalism could be used to calculate arbitrary higher mo- 
ments of the distribution of P e . 
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From the evolution equation 
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it is straightforward to show that 
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As before, we define the density matrix p, which now contains the square of 
the neutrino survival probability, and use the formalism of Section El with 
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(26) 



where the subscripts r and i denote the real and imaginary part of the quan- 
tity. 

We again solved the system numerically, assuming fluctuations of the 
same form as in Section |3|. The results are presented in Figures 7-9, for 
the same mixing parameters and values of (/3} rms as in Figure 1-3. The 
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quantity plotted is cr = yV, rather than the variance itself. For the smaller 
angles, there is a dramatic peak in a around Am 2 /E = 1.6 x 10~ 5 eV 2 /MeV, 
corresponding to neutrinos being produced in the resonance region (since 
we assume that all neutrinos are produced at the center of the sun). The 
results for sin 2 29 = 0.7 show a similar increase around Am 2 /E = 1.6 x 
10~ 6 eV 2 /MeV, which again corresponds to production inside the resonance 
region. It does not go back down, since the resonance region will include the 
core for any Am 2 /E larger than that. 

As the rms noise amplitude increases, the variance seems to saturate at 
a value of 1/12, independent of the mixing parameters. Two points are rele- 
vant to this phenomenon. First, a uniform probability distribution bounded 
between zero and one has the same mean and variance. Second, fluctuations 
will tend to suppress the oscillation of (P e ) as a function of time, since it 
becomes a superposition of oscillations with differing phases. Therefore, we 
may set the time derivatives on the left-hand sides of Eqs. ( |T8D and (|2B[) to 
zero. Assuming that k goes to zero smoothly, we are led to the trivial solu- 
tion r = for Eq. (|18D , and s = 1/6 for Eq. (|26|), where we have taken the 
limit k = 0. This yields a variance of 1/12 and a mean of 1/2, as expected. 

In order to illustrate the magnitude of the effect, in Figures 10-12, we plot 
the mean survival probability plus and minus the width a of the distribution. 
Note that, even though we did not calculate the higher moments, the mean 
survival probability plus and minus a is bounded between zero and one, 
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suggesting that the distribution of P e is not very skewed. 

5 Discussion 

The variance is a potentially important tool for the exploration of solar den- 
sity fluctuations on some time scales. Fluctuations on the time scale of a ra- 
diochemical experiment's run would broaden the distribution of count rates, 
since different runs would have different survival probabilities. None of the 
experiments currently operating has noted a broader distribution of rates 
than expected [13], [L4|, suggesting that the neutrino data will probably 
limit, rather than measure, such fluctuations. For the favored, small-angle 
solution, the variance is strongly peaked when neutrinos are produced near 
the resonance. This suggests that the finite radial distribution of neutrino 
production will have an important effect, likely extending the peak to lower 
values of Am 2 /E, since those correspond to resonance farther out in the 
sun. Further, it means that density fluctuations will affect pp neutrinos 
more strongly than other neutrinos, so that gallium experiments may put 
the strongest limit on fluctuations. 

To develop a very rough estimate of the limit on fluctuations, we note 
that their effect should become noticeable when the ratio a/(P e ) becomes 
comparable to the relative la experimental uncertainty. To simplify the 
argument, we assume that the signal observed at the gallium experiment 
consists only of ~ 0.3 MeV pp neutrinos. (Standard MSW analyses indicate 
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the near-complete suppression of other neutrinos.) As an example, we will 
consider the GALLEX results. GALLEX has a la uncertainty of approxi- 
mately 13% and an energy-averaged survival probability of about 60% |jl4f| . 
A 1% density fluctuation on the time scale of a GALLEX run, 20 to 28 days, 
with sin 2 29 = 0.01 and Am 2 /E = 1.63 x 10" 5 eV 2 /MeV (chosen to give a 
60% survival probability), has aj (P e ) = 0.15. It is likely, then, that a careful 
study could rule out fluctuations on that level. 



In the future, real-time, high-statistics detectors such as HELLAZ JL6 
and Borexino |I7] could be used to investigate fluctuations on shorter time 
scales. In particular, a helioseismological g-mode oscillation could leave a 
signature in the neutrino data. Unfortunately, the current high-statistics 
experiments SNO [18| and Super-Kamiokande [19| probably cannot be used 
to probe fluctuations in this way, as they are not sensitive to low-energy 
neutrinos. 

A proper averaging over neutrino production location would be extremely 
cpu-intensive. In this regard, approximate analytic techniques to compute 
the moments of the density matrix, such as that of Ref. |§, may be very 
useful in analyzing data. 
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Appendix 

Here we summarize the application of logarithmic perturbation theory to a 



two level system pOfl . Consider 



subject to the initial condition ipi(t = 0) = 1 and ^{t = 0) = 0. Defining 

-| <»> 

we find that the quantity z satisfies the Riccati equation, 

iz = -2Az + gC(l-z 2 ). (29) 

Expanding z in powers of g 

z = z + gz x + g 2 z 2 + ... , (30) 

inserting into Eq. (|2"9"D, and equating powers of g, we find that 

z = -ige+nfiW* tdt'Cit') e -*tfW)*" + (g 3 ) (31) 
Jo 
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is the solution with the given initial conditions. In order to calculate the 
survival probability, IV'iI 2 , we observe that 

1 



l^il 2 



(32) 



1 + |z| 2 ' 

where we have used \ipi\ 2 + 1^2 1 2 = 1- It is easy to show that Eq. ( p9f ) can 
be rewritten as 

zllog(l + H 2 ) =gC(z*-z). (33) 
Eqs. (§2]) and fl33|) then yield 

r-T 

(34) 



2 = exp 


2g f 




Jo 



2g [ dtC(t)lmz 
Jo 



Finally, substituting the approximate solution, Eq. (|3lD , into Eq. (|34D , we 
arrive at the desired answer: 

2 



|V>i | 2 = exp 



" 9 2 
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dt C(t) e +2i Io A ^ dt ' 



+ 0(g 4 



(35) 



It usually is possible to calculate the integral in Eq. (|35|) within the sta- 
tionary phase approximation. If the quantities A(t) and C(t) are monotoni- 
cally changing, there is only one stationary point t a and we can approximate 
the integral above as 

2 C 2 (t a ) 



exp 



-ng 



A'(t a ) 



0(g 4 



(36) 



where the prime denotes the derivative with respect to time. In most cases, 
especially when A(t) fluctuates, there could be more than one turning point. 
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If two of the turning points are close to each other it would be necessary to 
employ the Airy uniform approximation ^TL\. In this case one gets 



1-01 
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(37) 



j\A'(t a )\ y/\A'(t b )\ 

where t a and t b are the turning points, Ai(x) and Ai'(x) are the Airy function 
and its derivative, and 



c 



3 r*» 1 2 / 3 
A(t) dt 



(38) 



If the turning points are far apart, one can simplify Eq. (|3~?D using the asymp 
totic expressions of the Airy function: 



\i>i | = exp 



-2irg' 



\A'(t a )\ 
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sin" I — 



4 Jt a 



(39) 



where we have assumed that the functions C(t) and A'(t) are very slowly 
changing in this interval, i.e. C(t a ) = C{t b ) and A'(t a ) = A'(t b ). 
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Figure Captions 

Figure 1. Mean survival probability for a Bahcall-Pinsonneault density 
profile and sin 2 29 = 0.01. The probability is plotted for (/3) rms equal to 0, 
0.01, 0.02, 0.04, 0.08, and 0.5, with a constant correlation length r = 10 km. 

Figure 2. As Figure 1, but with sin 2 29 = 0.1. The vertical dotted line 
indicates where the constraint of Equation (|22|) begins to break down; there, 
sm29Am 2 /E = 0.1 x (3.95 x 10" 5 ) eV 2 /MeV. 

Figure 3. As Figure 1, but with sin 2 29 = 0.7. The vertical dotted line is as 
in Figure 2. 

Figure 4. Mean survival probability for the Bahcall-Pinsonneault sun com- 
pared to an exponential density profile with the same central and edge den- 
sity, with sin 2 29 = 0.01. (a)-(e) show (/3) rms equal to 0, 0.01, 0.02, 0.04, and 
0.08, respectively. The solid line is the Bahcall-Pinsonneault sun, and the 
dashed line is the exponential. 

Figure 5. As Figure 4, but with sin 2 29 = 0.1. The vertical dotted line is as 
in Figure 2. 

Figure 6. As Figure 4, but with sin 2 29 = 0.7. The vertical dotted line is as 
in Figure 2. 
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Figure 7. a for a Bahcall-Pinsonneault density profile, with the same values 
of sin 2 29 and ({3) vms as Figure 1. 

Figure 8. As Figure 7, but with sin 2 29 = 0.1. The vertical dotted line is as 
in Figure 2. 

Figure 9. As Figure 7, but with sin 2 29 = 0.7. The vertical dotted line is as 
in Figure 2. 

Figure 10. Mean survival probability plus and minus a for the Bahcall- 
Pinsonneault density profile with sin 2 29 = 0.01. (a)-(d) show ((3) rms equal 
to 0.01, 0.02, 0.04, and 0.08, respectively. 

Figure 11. As Figure 10, but with sin 2 29 = 0.1. 
Figure 12. As Figure 10, but with sin 2 29 = 0.7. 
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Figure 6e 
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